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ABSTRACT 

We investigate the environmental dependence of stellar population properties of galaxies in 
the local universe. Physical quantities related to the stellar content of galaxies are derived 
from a spectral synthesis method applied to a volume-limited sample containing more than 
60 thousand galaxies (0.04 < z < 0.075; M r ^ -19.9), extracted from the Data Release 4 
of the Sloan Digital Sky Survey (SDSS). Mean stellar ages, mean stellar metallicities and 
stellar masses are obtained from this method and used to characterise the stellar populations 
of galaxies. The environment is defined by the projected local galaxy density estimated from 
a nearest neighbour approach. We recover the star formation-density relation in terms of the 
mean light-weighted stellar age, which is strongly correlated with star formation parameters 
derived from Ha. We find that the age-density relation is distinct when we divide galaxies 
according to luminosity or stellar mass. The relation is remarkable for galaxies in all bins of 
luminosity. On the other hand, only for an intermediate stellar mass interval (associated to a 
transition in galaxy properties) the relation shows a change in galaxy properties with environ- 
ment. Such distinct behaviours are associated to the large stellar masses of galaxies with the 
same luminosity in high-density environments. In addition, the well known star formation- 
density relation results from the prevalence of massive systems in high-density environments, 
independently of galaxy luminosity, with the additional observed downsizing in galaxy for- 
mation, in which the star formation is shifted from massive galaxies at early times to low-mass 
galaxies as the universe evolves. Finally, our results support that a natural path for galaxy evo- 
lution proceeds via a nurture way, in the sense that galaxy evolution is accelerated in denser 
environments, that took place mainly at high-redshifts. 

Key words: galaxies: evolution — galaxies: formation — galaxies: fundamental parameters 
— galaxies: stellar content — stars: formation 



1 INTRODUCTION 

Galaxy environment has a fundamental role in the evolutionary 
paths followed by galaxies. It is well known that galaxy popu- 
lations change according to the environment in which they are 
found. The most classical evidence of such env ironmental de- 
pendence is the morphology-densi ty relation (e.g. |Pressleall980l : 
IWhitmore. Gilmore & Joneslll993r) . The fraction of galaxies with 
distinct morphological types (essentially spirals, lenticulars and el- 
lipticals) correlates strongly with the local galaxy density, with 
high-density environments being populated preponderantly by 
early-type galaxies. This is related to the dependence of the fraction 
of star-forming galaxies with the environment, based in the pres- 
ence or not of emission lines, such as Ha and [O I I ], in galaxy spec- 
tra (e.g. lCarter et al.l200ll ; lHashimoto et all 1 9981 ; iMateus & Sodrel 
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1 20041) . Such dependence is also intrinsically r elated to the reduced 
gaseous content of galaxies in denser regi ons dSolanes et alJI 19961 : 
iBravo-Alfaro et al.ll2000l : lGoto et al]2003l) . 

Recent studies have found that the star formation rate (SFR) 
of galaxies is the most sensitive parameter to galaxy environ- 
ment, de clining strongly in high-density regions of galaxy clus 
ters (e.g. lLewis et alj|2002t iGomez et al j|2003l : iRines et 



axy ctus - 
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whereas structural para meters are less affected by environment 
JKauffmann et al.ll2004) . The relation environment-SFR has also 
an additional de pendence on galaxy luminosity, being steeper for 
fainter galaxies JTanaka et alj|20041) . However, when one restricts 
the analysis to star-forming galaxies the median SFR of this class 
of objects seems to be unaffected by environment, although the 
fraction of such ga laxies decreases with increasing local density 
dBalogh et al]|2004h . 

The most plausible way to account for these trends is to 
assume that galaxy properties (mainly those related to star for- 
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mation and gas properties) are affected by environment through 
well known physical mechanisms acting on galaxies. This path, 
linked directly to the environment, gives origin to a nurture per- 
spective for galaxy evolution. In fact, several physical mechanisms 
were already proposed and studied to account for the evolution- 
ary trends discussed above. Interactions between the intragalactic 
and intergalactic mediu m, including gas removal and evaporation 
{ram pressure stripping; Gunn & Gotull972l : iFuiita & Nagashimal 
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1200 lh , and the suppression of the accretion of 



gas-rich materials in the neighbourhood of the galaxy (starvation ; 
lLarson, Tinslev & Caldweiilll980l : iBekki. Couch & ShiovdlJOOih . 
are typical mechanisms claimed to affect star formation properties 
of galaxies infalling onto denser regions. 

Another path by which galaxy evolution proceeds is related 
to the initial conditions established during the formation epoch of 
galaxies, which could, in principle, account for the relations be- 
tween galaxy properties and environment. This path gives origin to 
a nature perspective driving galaxy e volution. As expected from a 
biased galaxy formation scenario (e.g. lCen & Ostrikeiill993l) . mas- 
sive galaxies are formed earlier preferentially in high-density re- 
gions; in opposition, low mass galaxies would be formed later with 
a smooth distribution in the density field. Thus, in this scenario 
there is, naturally, an expected relation between age of formation 
(or mass) and environmental density. In fact, hierarchical galaxy 
formation m odels successfuly re produce the morphology-density 
relation (e.g. iBenson et alj|200ll) , demonstrating that it is present 
even if the physical mechanisms mentioned above are not consid- 
ered. 

Additionally , studies based on high-redshift galax ies (e.g. 
ICowie et al.ll996l : lKodama et al.l2004l:|juneau et alj2005l and oth- 
ers) and also on local galaxies (Kauffmann et al. 2003a; Heav- 
ens et al. 2004; Mateus et al. 2006, hereafter SEAGal II) have re- 
vealed that the existence of a 'downsizing' in galaxy formation is 
extremely important in analysis involving the star formation prop- 
erties of galaxies. These results suggest that massive galaxies have 
stopped to form stars at earlier times, with low mass systems com- 
prising a large fraction among galaxies with ongoing star forma- 
tion. These tren ds are also recovered by recent high-resolution 
simulations (e.g. IWeinberg et alj|2004j) and semi-analytic m odels 
for galaxy formation (e.g. lMenci et alj|2005i ; lde Lucia e t al. 2006). 
However, the physical origin of the downsizing is still a subject 
of debate (e.g.|de Lucia et al.l2006l ; IScannapieco. Silk. & Bouwensl 
l2005l : lBundvetalJl2006h . ~ 

Here we will shed some light into the discussion concerning 
galaxy formation and evolution by investigating the role of en- 
vironment on the stellar population properties of galaxies in the 
nearby universe. Other works have increas ed our understanding on 
this issue through distinct approaches (e.g . [Kauffmann et al. 2004; 
iThomas et al.ll2005l : |Poggianti et alj|2~005l : IClemen s et al. 2006). In 



this paper, we will use the results of the application of a spectral 
synthesis method to a volume limited sample of SDSS galaxies. 
This method is able to recover important physical properties of 
galaxies from their spectra, including mean stellar ages, mean stel- 
lar metallicities and stellar masses, among others. The details of 
such approach were discussed by Cid Fernandes et al. (2005) (here- 
after SEAGal I) and its results have been recently used in SEAGal 
II, where we have studied the bimodality of the galaxy population. 
In this work we show that the relations between galaxy prop- 
erties and environment, well characterised by the age-density re- 
lation, are distinct when we divide galaxies according to luminos- 
ity or stellar mass. This is a result of the different mass-to-light 
(M^ J L) ratios of galaxies in distinct environments. In dense re- 



gions galaxies of same luminosity tend to have higher M*/L ra- 
tios. We argue that a natural path emerges from these results but 
via a nurture way in the sense that massive galaxies are formed 
preferentially in high-density regions (where the evolution was ac- 
celerated) in the high-z universe (when the star formation rate was 
high), since their stellar populations also tend to be the oldest ones 
and they inhabit the densest regions we observe in the local uni- 
verse. 

The paper is organised as follows. Section 2 describes the data 
used in this work, the galaxy sample definition and a brief overview 
of the spectral synthesis method used in our analysis. Section 3 de- 
scribes some parameters related to star formation activity in galax- 
ies that will be discussed in this work, and Section 4 defines the 
environmental parameter, namely the projected local galaxy den- 
sity. In Section 5 we investigate the environmental dependence of 
the stellar population properties of galaxies. Finally, in Section 6 
we discuss the implications of our findings in the context of galaxy 
evolution, and in Section 7 we summarise the main results of this 
work. 



2 THE DATA 

The data used in this work were extracted from the Sloan Dig- 
ital Sky Survey available publ icly in its Data Release 4 (DR4; 
lAdelman-McCarthv et al.l2006l) . The SDSS is an ambitious project 
planned to study the large scale structure of the universe as well 
as other relevant subjects of extragalactic astronomy. The survey 
aims to cover one-quarter of the entire sky (mainly in the North- 
ern Hemisphere) obtaining photometric information in five optical 
bands (u, g, r, i, and z) for millions of objects and spectroscopic 
data for about 10 6 galaxies brighter than r — 17.77. 

In this section we describe the volume-limited sample of 
SDSS galaxies used in our analysis and the application of a spectral 
synthesis method to the spectroscopic data aiming to obtain some 
fundamental physical properties of galaxies in the local universe. 
In this work we use the following values of the cosmological pa- 
rameters: Ho = 70 km s _1 Mpc -1 , Q M = 0.3 and fi A = 0.7. 



2.1 Galaxy sample 

We built a volume limited sample of galaxies from the SDSS DR4 . 
This sample follows the same criteria used bv lStrauss et alJJ2002l) 
to define the Main Galaxy Sample, consisting of galaxies with r- 
band Petrosian magnitudes r ^ 17.77 and r-band Petrosian half- 
light surface brightnesses /xso $J 24.5 mag arcsec 2 . We have se- 
lected galaxies with redshifts in the range 0.04 < z < 0.075 with 
r-band extinction-corrected absolute magnitudes M r ^ —19.9, 
corresponding roughly to M* + 1.5. The absolute magnitudes 
adopted here are fc-corrected by using the kcorrect v3.2 code pro- 
vided bv lBlanton et alj J2003bh . The lower redshift limit was cho- 
sen here to avoid the undesirable presence of aperture effects in our 
analysis which are related to the use of small 3" fibre s in the SDSS 
to obtain the spectroscopic data (Kewlev et al.ll2005l) . We have de- 
tected a minor fraction of galaxies with multiple spectral data, from 
which we have considered only those with best signal-to-noise ra- 
tio (S/N) in the (/-band. These selection criteria result in a volume 
limited sample containing 63659 galaxies. 
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Figure 1. Mass-luminosity relation for (a) all galaxies in the 20k sample, and for (b) late-type galaxies with D n (4000) < 1.67 and (c) early-type galaxies 
with D n (4000) > 1.67. Solid lines are the median values (thick lines) and 10-90 percentiles (thin lines) for the relations. The luminosity limit for our sample 
(log L r / Lq = 9.80) is shown as vertical dashed lines, while the mass values above which our sample will be considered complete are shown as horizontal 
dashed lines. The numbers at the top-right are the count of galaxies in each panel and the Spearman correlation coefficient. 



2.2 Overview of the spectral synthesis method 



We used th e starligh t code (described in detail in lSEAGall and 
updated in ISEAGal DT) to obtain physical parameters of galaxies 
directly from their spectra (essentially from continuum and absorp- 
tion lines), starlight fits an observed spectrum with a combina- 
tion of N* Simple Stel lar Populations (SSP) from the evolution- 
ary synthesis models of iBruzual & Charlotl ( 120031) . Here, we used 
a base containing N* = 150 elements, spanning 6 metallicities - 
Z = 0.005, 0.02, 0.2, 0.4, 1 and 2.5 Z e - and 25 ages, between 
1 Myr and 18 Gyr. Stellar extinction is modelled as due to fore- 
round dust, with the extinction law o flCardelli. Clayton. & Mathij 
19891) with Rv = 3.1, and parametrised by the V-band extinction 
Ay- Line of sight stellar motions are modelled by a Gaussian dis- 
tribution centred at velocity «* and with dispersion a*. Regions 
arou nd emission lines and bad pixels are excluded from the anal- 
ysis. ISEAGallJ describes numerical and other technical aspects of 
the code. 



We apply our synthesis method to the galaxy sample described 
in the last section. The SDSS spectra cover a wavelength range 
of 3800-9200 A, have mean spectral resolution A/AA ~ 1800, 
and were taken with 3 arcsec diameter fibres. Prior to the spectral 
fits, the s pectra are corrected for Galacti c extin ction with the maps 
given by Schlegel, Finkbei ner. & Davis! J1998I) and the extinction 



law of ICardelli et al.l dl989l ), shifted to the rest-frame and resam- 



pled from 3400 to 8900 A in steps of 1 A. 

The result of the procedure is a list of parameters for each 
galaxy. The most important for the discussion in this paper are: 
a) the dust-corrected stellar mass inside the fibre; the total stellar 
mass (M+) is computed a posteriori after correcting for the fraction 
of luminosity outside the fibre assuming that the galaxy mass-to- 
light ratio does not depend of the radius; b) the mean stellar age, 
weighted either by light, (logt*)_L, or by mass {log£*)jvr; c) the 
mean stellar metallicity, light and mass-weighted, \og{Z- k )L and 
log{Z*)M, respectively; d) the V-band stellar extinction A v ; ande) 
the velocity dispersion, <r*. Besides, by subtracting the synthesised 
spectrum from the observed one, we get a "pure-emission" residual 
spectrum, which is useful for analysis of the galaxy emission lines. 



2.3 Completeness issues 

The volume-limited sample adopted in this work secures a high 
completeness for galaxies brighter than M r ^C —19.9. However, 
in this work we are also interested in the properties of galaxies di- 
vided according to stellar mass, thus a high completeness in this 
parameter is needed in order to avoid any bias originated from the 
sample selection procedure. 

We investigate this issue with the help of the mass-luminosity 
relation of galaxies from a flux-limited sample containing 20000 
objects (hereafter referred as '20k sample') selected at random. The 
total stellar masses are computed after correcting the stellar masses 
obtained from the spectral synthesis by aperture effects. As for the 
main sample, this is done by dividing the stellar masses by (1 — /), 
where / is the fraction of the total galaxy luminosity in a given 
band outside the fibre. Here this correction is made by using total 
and fibre Petrosian magnitudes in the z-band. 

Fig. Q] shows the mass-luminosity relation for all galaxies in 
the 20k sample, and for early and late-type galaxies distinguished 
according to the 4000 A break index measured following the def- 
inition given by iBalogh et al.l d 1 9991) and referred as D n (4000). 
In SEAGal II we have shown that the optimal value to separate 
early and late galaxy types is D n (4000) = 1.67. In this figure we 
also show the median values (solid lines) and the 10-90 percentiles 
(dashed lines) of the stellar mass in bins of galaxy luminosity. 

In order to select a subsample of galaxies with high complete- 
ness in stellar mass we need to compute the value of M* above 
which it will be complete. Adopting the galaxy luminosity limit for 
our sample (M r ^ —19.9, corresponding to log L t /Lq ^ 9.8, 
shown in Fig. Q] as vertical dashed lines), and considering the up- 
per 90 percentile of the distributions for each galaxy type shown in 
Fig.Q] we found that our galaxy sample will be nearly complete for 
late-type galaxies with Af* ^ 1.8 x 10 M©, and for early-type 
galaxies with Af* ^ 2.8 x 10 10 M©. By using these limits (shown 
as horizontal dashed lines in Fig. [TJ, we then selected a subsam- 
ple containing 49453 objects, with a high completeness in stellar 
mass. Thus, when studying relations involving M*, we will adopt 
this subsample in order to avoid any bias introduced by our sample 
selection procedure. It is convenient to summarise here the sam- 
ples that will be used in the rest of the paper. In Table Q] we show 
the basic statistics of the volume-limited and stellar mass-limited 
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Figure 2. Relations among the tracers of star formation activity based on the Ha emission line. Median values of EW(Hq) as a function of SFR(Ho) (left 
panel) and SSFR(Hq) (right panel) are shown in bins of r-band galaxy luminosity chosen to contain the same number of objects. In this and in the following 
figures, the median values for each luminosity bin are shown as different lines. The legend in the top-left corner shows the range of each bin in which the 
median values have been evaluated. The numbers at the top-right are the count of galaxies in each panel and the Spearman correlation coefficient. 
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Figure 3. Median values of EW(Ha) and SSFR(Hq) as a function of light-weighted mean stellar age in bins of galaxy luminosity. The numbers at the top-right 
are the count of galaxies in each panel and the Spearman correlation coefficient. 



Sample 



Number Per cent 



Volume-limited 
Mi, -limited 
Star-forming galaxies 



63659 100.0 

49453 77.7 

26540 41.7 



Table 1. Summary of the samples adopted in this work. 



3 THE STAR FORMATION ACTIVITY IN GALAXIES 

In this section we describe some usual tracers of star formation 
activity in galaxies based on the Ha emission line, and how they 
are related to spectral synthesis products. Our main intention here 
is to define these quantities and show that the light-weighted mean 
stellar age behaves as a good tracer of current star formation in 
galaxies. 



samples, and also of a subsample of star-forming galaxies that will 
be analysed in the next section. 



3.1 Emission line related products 

The equivalent width (EW) of the Ha emission line is an useful 
parameter when studying the star formation properties of galax- 
ies. This quantity is related to the ratio between the amount of star 
formation occurring in a galaxy in the last ~ 10 7 yr and its past in- 
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tegrated star formation history. Since the measurement of EW(Ha) 
does not require flux-calibrated spectra nor intrinsic reddening cor- 
rection, it has been widely used as an indicator of recent star forma- 
tion activity in local galaxies, for instance in works related to the 
envi ronmental dependence of star formation properties of ga laxies 
(e.g. lGomez et al.ll2003l : Lewis et allioollRines et alj20O5t) . 

A more direct and usefu l tracer of star fo rmation rate in galax- 
ies is the Ha luminosity (e.g. lKennicuti 19980 . When appropriately 
corrected for underlying absorption and dust extinction, and in the 
case of fibre-based spectra also for aperture effects, the Ha SFR es- 
timate is f ound to be co nsistent with radio, far-infrared, and u-band 
estimates JHopkins et a l. 2003). In our case, the underlying stellar 
absorption is automatically subtracted when computing the resid- 
ual spectra (as discussed in Section l2~2l . The Ha SFR, corrected 
by aperture e ffects, is then sim ply derived by using the expression 
adapted from lKennicuttl i 19981) for a Chabrier IMF (0.1 -100 M Q ) 
and with the prescriptions given bv lHopkins et al.1 ( 120031) : 



SFR Ha (M yr _1 ) = 5.22 x 10~ 42 L(Ha) 10 



-0.4(r-7- fibro ) 



(1) 



where L(Ha) is the observed luminosity of Ha (in ergs s _1 ) cor- 
rected by nebular extinction with the Ha/H/3 intrinsic intensity 
ratio, r is the r-band Petrosian magnitude representing the total 
galaxy flux, and rfibrc is the r-band fibre magnitude^ Note that 
since we have corrected the Ha luminosities for extinction with the 
help of the Ha/H/3 ratio, only galaxies with these two lines mea- 
sured with significant S/N ratio will have their SFR(Ha) com- 
puted. Consequently, only galaxies with higher levels of star for- 
mation activity will have their star formation rates evaluated. 

Another interesting quantity related to star formation in galax- 
ies is the SFR per unit stellar mass, or specific star formation rate 
(SSFR(Ha)), which is associated to the strength of the current burst 
of st ar formation relative to the underlying stellar mass of a galaxy 
(e.g. lGuzman et al.ll 19971) . In our case, this parameter is easily ob- 
tained by dividing the SFR(Ha) (in M© yr _ 1 ) by the dust-corrected 
stellar mass, M*. 

In order to illustrate the quantities described above, in Fig. [2] 
we show the relation between the EW(Ha) and the SFR(Ha), and 
that between the EW(Ha) and the specific Ha SFR. This figure 
also shows the median values of each relation in bins of r-band 
galaxy luminosity. It is important to stress here that we only con- 
sider galaxies wi th Ha mea sured in emission with a S/N greater 
than 3 (following lSEAGal J) . In addition, we also restrict this analy- 
sis to galaxies with Ha emission coming from normal star-forming 
regions by exc luding host s of active galactic nuclei in the same 
way as done bv lSEAGal 111 In Fig. [2] we note a clear dependence on 
galaxy luminosity for the relation between EW(Ha) and SFR(Ha), 
in the sense that for a given value of EW(Ha), brighter galaxies 
have higher SFRs. This luminosity dependence disappears when 
dealing with the specific SFR. 



3.2 Light-weighted mean stellar age as an indicator of star 
formation 

The mean light-weighted stellar age of a galaxy is related to the 
formation epoch of massive and bright stars, frequently associated 



The latter term in this equation is used to correct the SFR estimate by 
aperture effects by assuming that the Ha luminosity inside the fibre is char- 
acteristic of the whole galaxy, and that the continuum flux at the Ha wave- 
length is well represented by the flux at the e ffective wavelength of the 
r-band filter (see details in lHopkins et aU2 003). 




Figure 4. Comparison between local and global definition of environment. 
The local galaxy density (£io) is plotted as a function of the cluster-centric 
radius (r/r v ; r ; a i). Median values of the relation are shown as a solid line, 
and their respective quartiles as dashed lines. In the top-right corner is 
shown the number of galaxies and the Spearman correlation coefficient for 
the relation. 



to starbursts. Hence, low values of (logt*)i reflect the recent ac- 
tivity of intense star formation operating in a given galaxy, whereas 
high values of it indicate that most of the galaxy light comes from 
older stellar populations, with less significant, or even none, recent 
star formation episodes. As this quantity is obtained by our spec- 
tral synthesis approach in a robust way we will inspect its ability in 
tracing star formation activity of galaxies. 

In Fig.[3]we show the relation between (log£*)z, and the pa- 
rameters discussed previously: EW(Ha) and SSFR(Ha). We do not 
show the relation between (log /*)l and SFR(Ha) since this latter 
parameter has a dependence on galaxy luminosity, related to the 
amount of stellar mass recently formed in a galaxy. As in Fig. [2] 
these relations are also shown in bins of galaxy luminosity (i.e., it is 
an extensive quantity, whereas (log/*)z,, EW(Ha) and SSFR(Ha) 
are all intensive quantities). We clearly note that this mean stel- 
lar age is very sensitive to star formation activity estimated via 
Ha products. The relations with EW(Ha) and SSFR(Ha) show 
high Spearman correlation coefficients (—0.58 for both relations). 
Hence, the (log t *}l figures out as a good tracer of recent star for- 
mation activity in galaxies, and since we have measured it for all 
galaxies in our sample (unlike Ha based indicators) it is well suited 
for analysis involving the whole galaxy population, even for galax- 
ies without emission lines. 



4 THE ENVIRONMENT 

In this work we use a conventional approach to dealing with galaxy 
environment. We adopt a non-parametric method to determine the 
local number density of galaxies, based on the fcth nearest neigh- 
bour density estimator (fcNN; see, for instance, iFukunagal i 19901) 
for a statistical description of such approach). This method fixes a 
value for k and lets the volume V(r), centred on a given object and 
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extending to its fcth nearest neighbour, be a random variable. This 
volume is large in low density regions and small in high density re- 
gions. In this way, this method provides a spatial description of the 
densit y field, which has been recently explored bv lMateus & Sodrel 
d2004h in a study of field galaxies. 

When one goes to denser environments, velocity dispersion is 
high and neglecting peculiar velocities in distance estimates could 
underrate local densities determined with the fcNN method. In or- 
der to avoid this effect, a better way to derive local galaxy den- 
sity via nearest neighbours is by considering a projected distribu- 
tion of galaxies, instead of a spatial one. This procedure has been 
intensively used in st udies devoted to galaxy environments, since 
the classical wor k bvlDresslei dl980f) to more recent works based 



20021) and SDSS (e.g. iGomez et al.l 



on 2d FGRS (e.g. iLewis et al. 
120031) data. 

A more common approach when dealing with large red- 
shift survey data is to define a limit in redshift space around a 
given galaxy, for instance in th e range ±1000 km s -1 , to find 
its fcth nearest neighbou r (e.g. iBlanton et al.l l2003al : iGoto et al.l 
l2003l ; lBalogh et al.ll2004f) . This limit emulates a three-dimensional 
description of the density field by excluding galaxies in the 
foreground/backgro und with an objective criterion. As noted by 
ICooper et alj fl2005t) . projected density estimates with limits in the 
line-of-sight velocities between ±1000 and ±1500 km s _1 are best 
suited for the description of a broad range of environments. 

Here we use the expression E^ = k/wr%, to estimate the local 
galaxy densitjQ where r^ is the projected distance to the fcth near- 
est neighbour of a given galaxy. We adopt the usual limit in redshift 
space of ±1000 km s _1 . Additionally, in our density estimates, we 
have prevented an incorrect determination of £& due to border ef- 
fects by excluding galaxies whose fcth neighbours have projected 
distances greater than the distance of the galaxy to the closest bor- 
der of the survey region or of the sample volume. Note that in the 
search of neighbours and in the density estimates we have used all 
galaxies in the volume limited sample independently of the quality 
of their spectroscopic data. 

We have computed our local galaxy density estimates with the 
values fc = 5 and fc = 10 in order to assess the influence of the 
choice of this parameter in the results presented in this work. How- 
ever, we have found that none of the results discussed in the next 
sections depend on the values of the fc tested here. In the rest of this 
paper we will adopt the local galaxy density as evaluated from the 
distance to the tenth neighbour, £10. 

We have compared our local density estimates defined by £10 
to a global parameter: the distance to the nearest cluster centre 
from each galaxy in our sample. We searched for galaxy clusters 

S iresent in our sam ple with the help of the C4 Cluster Catalogue 
Miller et alJl2005D . which contains 748 clusters identified in the 
spectroscopic sample of the SDSS Data Release 2 (DR2). This cat- 
alogue is about 90 per cent complete and 95 per cent pure above 
M200 = 1 x 10 14 IT 1 M and within 0.03 < z < 0.12, thus cov- 
ering the redshift range of our sample. Therefore, for each galaxy in 
our sample we can compute its distance to the nearest cluster cen- 
tre, defining in this way a cluster-centric radius normalised by the 
cluster virial radius. In Fig. [4] we show the relation between our 
density parameter (£10) and the cluster-centric radius (r/r v i r iai) 
only for galaxies located in the regions covered by the DR2 data. 



2 Note that in statistics, fc — 1 should replace fc in the numerator of the 
expression for S ^ in order to achieve an unbiased estimate of local density 
JFukunagall990l) . 



There is a clear correlation between the two quantities, mainly for 
r/r v irial i$ 5, above which we note that the median values of local 
galaxy density are almost constant. 

A possible source of bias in the estimation of local galaxy den- 
sity in the SDSS data is due to a technical limitation in the proce- 
dure of assignment of fibres to spectroscopic targets closer than 
55" in the sky. This problem, referred as 'fibre collision', becomes 
very important in denser environments and about 8-9 per cent of 
select ed objects are unobserved for this reason dStoughton et al.l 
l2002h . 

In order to assess the significance of this bias in our results 
we have searched in the DR4 photometric catalogue for galaxies 
closer than 55" of each galaxy in our spectroscopic selected sam- 
ple, which actually have been rejected in the fibre assignment pro- 
cedure. Note that these missing galaxies only affect our local galaxy 
density estimate if they have radial velocities close to those of the 
selected galaxies. However, as they do not have any redshift infor- 
mation we have overestimated the effects of the fibre collision bias 
by supposing that all such missing 'neighbours' are at the same 
redshift as the selected galaxies. From the total of 63659 selected 
galaxies in our volume limited sample, we have found about 12 
per cent of galaxies with missing 'neighbours' closer than 55". We 
have then evaluated our density estimates by considering that these 
missing objects were in fact close neighbours. After redoing our 
analysis with these new densities we concluded that none of the 
results shown in this work is affected by the bias due to fibre colli- 



5 THE ENVIRONMENTAL DEPENDENCE OF GALAXY 
PROPERTIES 

In this section we investigate how some physical properties related 
to the stellar populations of nearby galaxies behave as a function 
of the local galaxy density. We begin by investigating the depen- 
dence of the fraction of galaxies of distinct mean stellar ages and 
masses with the environment, which can be linked to the classi- 
cal morphology-density relation. The environmental dependence 
of star-formation activity in galaxies, specially if estimated through 
the light-weighted mean stellar age parameter, is further investi- 
gated. We also analyse how the metallicities of galaxies vary along 
the density range studied in this work. Finally, the role of stellar 
mass and luminosity in driving the star-formation-density relation 
is investigated here, mainly through a principal component analysis 
applied to our parameter set. 



5.1 Galaxy fractions 

It is well established that the fraction of galaxies classified accord- 
ing to distinct morphological types is strongly related with the en- 
vironment, with high-density regions typical of galaxy clusters be- 
ing populated essentially by elliptical galaxies, i n contras t to th e 
high fraction of spirals found in the field (e.g. iDresslerl ll980T) . 
This galaxy type segregation is investigated here in terms of two 
physical galaxy properties: the mean light-weighted stellar age and 
the stellar mass. SEAGal II studied these parameters in light of 
the bimodality observed in galaxy populations and found that the 
mean stellar age of galaxies is the primary responsible for the di- 
chot omy between blue and red galaxies seen in the local universe 
(e.g. lStrateva et al.ll200ll : lBaldrv et alJl2004l) .Thev also found that 
the stellar mass has an additional role in the sense that, blue, star 
forming galaxies are preferentially of low mass. 
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Figure 5. Left: Fraction of galaxies in bins of light-weighted mean stellar age (each one containing the same number of galaxies) as a function of the local 
galaxy density. Right: The same, but now for galaxies in bins of dust-corrected stellar mass. The numbers at the top-right are the count of galaxies in each 
panel. 



In Fig.[5]we show how the fraction of galaxies in bins of light- 
weighted mean stellar age (left panels) and stellar mass (right pan- 
els) correlates with the local galaxy density. Galaxies are grouped 
in four bins of (log £*}l and log M*/Mq, each one containing the 
same number of objects. It is clear from Fig.[5]that high-density en- 
vironments are populated by galaxies with older stellar populations 
and larger stellar masses, with the fraction of these objects increas- 
ing as the environment becomes denser. The transition in galaxy 
fractions occurs at a local density of S* ~ 0.7 h?o Mpc -2 , above 
which the population of massive galaxies with older stellar popu- 
lations begins to dominate. This transition density corresponds, on 
average, to about 2-3 virial radii from cluster centres, as can be de- 
duced from Fig. [4] Additionally, we note that in the case of A/*, the 
fractions only change for the two extreme bins, whereas intermedi- 
ate ones possess constant fractions along £10 . It is worth stressing 
that the transition shown in Fig.[5]is not th e same as the "brea k" at 
a characterist i c gala xy density reported bv lLewis et al.l ( 120021) and 
I Gomez et alj ( 120031) that we discuss below, since here we are in- 
vestigating only galaxy fractions as a function of the local galaxy 
density, instead of the SFR-density relation itself. 



5.2 Star formation activity 

Here we revisit the distribution of the Ha equivalent width as a 
function of the local galaxy density. In recent years, many works 
have been d e voted to the study of this relation. For instance, 
I Gomez et alj ( 120031) have investigated the environmental depen- 
dence of star formation rate in galaxies based on the Ha emis- 
sion line. Their main result was the finding of a break at a char- 
acteristic local galaxy density resulting in an abrupt decrease of 
the EW(Hq) (or S FR) of galaxies in denser regions. Recently, 
iTanaka et alj J2004D have shown that this break indeed occurs only 
for fainter galaxies (M* + 1 < M r < M* + 2), with brighter 
(M r < M* + 1) ones showing EW(Ha) only monotonically de- 
creasing with density. Here we use a more detailed approach to 
analyse this relation. 

In Fig. [6] we present the median values of EW(Ha) as a func- 
tion of local galaxy density. The median values of EW(Ha) are 



taken for galaxies grouped in bins of luminosity, log L v /Lq (left 
panel) and stellar mass, log M^/Mq (right panel). In this figure 
we show all galaxies with Ha emission line measured with S/N 
greater than 3, avoiding to include active nuclei hosts (c.f. Sec- 
tionO. This procedure is very robust since the line measurenments 
we have done make p ossible a detailed galaxy classification, as dis- 
cussed in lSEAGal 11 Moreover, with the help of spectral synthesis 
we can measure only the emission contribution of Ha line, avoiding 
to include in our analysis galaxies that would appear with negative 
values in distributions of EW(Ha ) or SFR(Ha), as s hown in other 
works (see, for instance, Fig. 4 of lGomez et alj2003f) . 

We note in Fig. [6] that the median values of EW(Ha) for 
galaxies in all luminosity bins decrease with galaxy density with- 
out showing any break (even if we conside r all quartiles of the 
distribution), confirming the result found by ITanaka et alj 120041) 
for brighter galaxies (our sample contains galaxies with M r < 
M* + 1.5). In this sense, the relation between star formation ac- 
tivity and local galaxy density for bright galaxies is independent 
of galaxy luminosity, since it exists in all bins of log L v /Lq . This 
trend is also seen in the right panel of Fig.[6] where galaxies are now 
divided according to their stellar masses. The absence of a break in 
Fig.[6]indicates that the relation between star formation activity and 
local galaxy density decreases monotonically over a w ide range of 
environments, as suggested bv lMateus & Sodrej J2004I) . 

We also investigate this issue by analysing the environmen- 
tal dependence of the star formation activity of galaxies spectrally 
classified a s star-forming, following the same definition as used in 
ISEAGal 11 These galaxies have low stellar masses (about 83 per 
cent of them have log M+/Mq < 10.67) and young mean stel- 
lar ages (more than 95 per cent of them have (log£*)x, < 9.53). 
In Fig. [7] we show the specific SFR (SSFR(Ha)) as a function of 
the local galaxy density only for the spectral class of star-forming 
galaxies. As in Fig. [6] galaxies are divided in bins of luminosity 
(left panel) and stellar mass (right panel). We note in Fig.|7]that the 
median values of SSFR(Ha) for star-forming galaxies have almost 
constant values along the density bins covered by our sample, im- 
plying that the specific star formation rate of star-forming galaxies 
does not depend on the local galaxy density, confirming the results 
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Figure 6. Median values of the Ha equivalent width as a function of local galaxy density, in bins of galaxy luminosity (left panel) and stellar mass (right 
panel). The numbers at the top-right are the count of galaxies in each panel. 



-9.9 






'-10.0 



13-10.1 
X 

t^-10.2 
CD 



-10.3 



Star-forming 

logL r 

9.80- 

9.90- 

- 10.02 

10.19 



N = 15228 



9.90 
10.02 
10.19 
10.98 




Star-forming 

log M, 

10.25 - 10.38 

- 10.38-10.51 

10.51 - 10.71 
10.71-11.46 



N = 8157 




^ \y v 



x b'ffll> W 



0.1 



1.0 

,2 A/ r -2l 



10.0 



0.1 



£ 10 [h /0 Mpc- 2 ] 



1.0 

.2 A/ r -2l 



10.0 



£ 10 [h /0 Mpc-1 



Figure 7. Median values of the specific SFR related to the Ha emission luminosity as a function of local galaxy density, in bins of galaxy luminosity (left 
panel) and stellar mass (right panel), only for galaxies classified as star-forming. The numbers at the top-right are the count of galaxies in each panel. 



found by other investigators (e.g. iBalogh et alj|2004l ; iRines et al.l 
120051) . 



5.2.1 Mean stellar ages 

An important issue concerning to the use of star formation trac- 
ers based on emission line measurements - like those related to 
the Ha emission line - is the absence of data for galaxies with lit- 
tle (or even none) recent star formation activity, since the emission 
lines (if they exist) in the spectra of these systems are too weak to 
be measured with appreciable signal-to-noise ratio. Additionally, 
in galaxies where a starburst is obscured by the presence of huge 
amounts of dust, the emission lines will also become very weak to 
be detected. Thus, results from environmental studies based solely 
on the EW(H«) or SFR(Ha) parameters could suffer from this kind 
of bias. In Section [3] we have shown that the mean light- weighted 
stellar age obtained by spectral synthesis has a good correlation 



with both EW(Hq) and SSFR(Ha), implying that this parameter is 
well related to the star formation history of galaxies. In fact, low 
values of (log £*)_£, reflect an intense activity of star formation oc- 
curring in a galaxy, whereas high values of it are associated to most 
of the galaxy light coming from older stellar populations. As we 
measured this parameter for all galaxies in our sample, it is inter- 
esting to investigate its dependence on the environment defined by 
the local galaxy density. 

In Fig. [8] we show the mean stellar age weighted by light as a 
function of the local galaxy density. The median values of (log £* ) l 
are taken for galaxies grouped in bins of luminosity and stellar 
mass. The behaviour of the mean ages of the stellar populations is 
distinct when we split galaxies according to L r or M* . In order to 
clarify these trends in terms of the bimodal character of the galaxy 
populations, we also show in this figure, as horizontal dashed lines, 
the value of (log^)^ which better separates early and late galaxy 
types (see lSEAGailil for details). 
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We note in the left panel of Fig. [8] that the median values 
of (logi*)i increase significantly as the environment becomes 
denser, independently of the galaxy luminosity, with the relation 
for fainter galaxies being steeper than that of brighter ones. The 
properties of the population of fainter galaxies also tend to change 
at Eio < 5 hro Mpc~ 2 , where the median values of (log £*)l cross 
the dividing line which distinguishes galaxies dominated by young 
stellar populations (late-types) and galaxies with older stars (early- 
types). On the other hand, the mean stellar ages of galaxies with 
different stellar masses (seen in the right panel of Fig. [8} behave in 
a distinct way with respect to the local galaxy density. Low-mass 
galaxies have median values of mean stellar ages slightly varying 
along the bins of density, showing low values of (log £*)i - char- 
acteristic of late-type galaxies - in all environments. However, this 
behaviour for the lowest stellar mass bin may be affected by incom- 
pleteness due to the low fraction of early-type galaxies selected in 
the range 10.25 < log M* < 10.48. The most massive galax- 
ies also show a slight change in the median values of (log£*)i, as 
galaxy density increases, but their values are higher everywhere. 
The steeper relation is seen for galaxies with intermediate stellar 
masses (~ 3 - 6 x 10 10 Mq), which shows a transition from late 
to early galaxy types at Eio ~ 0.3 — 1.0 h^o Mpc -2 . 

Therefore, an age-density relation is clearly obtained when 
galaxies are divided according to their luminosities. However, the 
relation is not so evident when they are grouped in bins of stellar 
mass. This indicates that the galaxy luminosity and stellar mass 
could be playing distinct roles in driving the age-density relation 
shown in this work. We will appropriately address this question in 
section[5/fl 



5.3 Stellar metallicities 

We now investigate the environmental dependence of another im- 
portant physical parameter of galaxies, the mean stellar metallicity. 
In Fig. [9] we plot log(Z*)_L and log(Z*) a/ versus local galaxy den- 
sity, for galaxies in bins of luminosity (left panels) and stellar mass 
(right panels). The relation between log(Z,)i and local density be- 
comes less steep as galaxy luminosity increases. The result of this 



trend is that fainter galaxies in low-density regions have subsolar 
metallicities, whereas in denser regions they have high metallicities 
similar to that of the brightest galaxies in our sample. For galaxies 
divided according to stellar mass, the trends are similar. The stellar 
metallicity of low mass galaxies increases as the environment be- 
comes denser, while massive galaxies have almost constant values 
of \og(Zi,)L along density bins. In the case of the mass weighted 
mean stellar metallicity (log(Z*)j\/) versus local galaxy density, 
the trends are similar, although less steep than those for log(.Z*)L. 
The main difference is that luminous (massive) galaxies also tend 
to have more metallic stellar populations in denser environments, 
as shown in the bottom panels of Fig. [9] 

It is interesting to note that the values of log(Z*)z, are re- 
lated to mean metallicities of the stellar populations which con- 
tribute significantly to galaxy light. Thus, one could expect that 
they are directly associated to recent episodes of star formation in 
galaxies. In this way, low-mass (or fainter) galaxies in denser envi- 
ronments, which are currently forming stars (following the results 
from Fig. [8), tend to have more metal-rich stellar populations than 
their counterparts in regions of lower densities. Low mass galax- 
ies have substantial amounts of gas to keep their high star for- 
mation rates. This gas could have been pre-enriched earlier pref- 
erentially in high-density environm ents, where the star formation 
activity was higher in the past (e.g. Madau. Ferrara & ReejuOOlc 



IScannapieco. Ferrara & Madaul I2002T) . Thus, one explanation is 
that these galaxies have already been formed in that denser environ- 
ments, but have not stopped to form stars until the present epoch. 
Thus, this result is additionally s upported by the observed 'down- 
sizing' in galaxy formati on (e.g. lCowie et al.ll 19961 ; iKodama et al.l 
120041 ; Uuneau et alj|2005l among others), where massive galaxies 
have formed most of their stellar masses at early times. However, 
the processes that have kept the continuous star formation in those 
low-mass galaxies are not well understood. 



5.4 The role of stellar mass and luminosity 

The relation between the mean stellar age of galaxies and their en- 
vironment (defined by the local galaxy density) shows that galaxy 
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Figure 9. Median values of light-weighted stellar metallicity as a function of local galaxy density for galaxies in bins of luminosity (left panel) and stellar 
mass (right panel). The numbers at the top-right are the count of galaxies in each panel. 



luminosity and stellar mass play distinct roles in defining the en- 
vironmental dependence of star formation properties of galaxies. 
In general, that dependence exists independently of L and only for 
an intermediate stellar mass range characterised by a transition in 
galaxy properties. 

Another way to visualise these trends is shown in Figs.llOland 
1111 where we show the distributions of (logt^)L for galaxies in 
vertical bins of local galaxy density, and horizontal bins of stellar 
mass (Fig. UOt and luminosity (Fig. lilt. Each bin in stellar mass 
(or luminosity) has been built to have the same number of objects, 
whereas the bins in log Eio are equally spaced. From the compar- 
ison of these figures, we note that the age distributions are distinct 
when galaxies are divided according to stellar mass or luminos- 
ity, mainly for objects in the extreme bins. In denser environments 
(corresponding to 0.91 < log Eio < 1.86), less luminous galaxies 
have a (log i*)_L distribution peaked at older ages, whereas in low- 
density regions the distribution is double-peaked, unveiling the bi- 
modal character of the galaxy population. On the other hand, in the 
case of galaxies divided according to stellar mass, the (log t* ) l dis- 
tribution for low-mass galaxies in denser environments is broader, 
without showing any peak. As local density decreases, the distri- 



butions are peaked at younger ages, indicative of the dominance 
of low-mass star-forming galaxies in low-density environments. At 
the other extreme, represented by very luminous and massive galax- 
ies, the distributions are similar, with the only difference being due 
to a larger fraction of galaxies with low values of (logi*)_L in the 
distributions for luminous galaxies. It is also interesting to note the 
differences in the distributions for two extreme environments. In 
denser regions, galaxies divided by luminosity show a conspicuous 
peak centred at older ages in all bins considered in Fig.QT] On the 
other hand, for galaxies divided by stellar mass there is a growth of 
the peak at (logi*)_t > 9.5 from the low-mass to the most mas- 
sive bin shown in Fig.QO] However, note that for the lowest stellar 
mass bin the absence of older early-type galaxies may be a reflex 
of our sample selection procedure to overcome the stellar mass in- 
completeness in our sample (Sectionf23). 

We also investigate the environmental dependence of a combi- 
nation of stellar mass and luminosity parameters, namely the mass- 
to-light ratio of galaxies in our sample as inferred by spectral syn- 
thesis. In Fig.[l2] we show the relation between the galaxy mass-to- 
light ratio (in r-band) and local density in different bins of r-band 
galaxy luminosity for (a) all galaxies in our sample, (b) early-type 
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galaxies with D n (4000) > 1.67 and (c) late-type galaxies with 
D n (4000) < 1.67. We note that, at fixed luminosity, galaxies in- 
habiting high-density environments tend to be more massive than 
their counterparts in regions of low density; the M*/L r ratio of 
these objects increases by 40 — 50 per cent from low to high-density 
regions. The relation is also steeper for the faintest bin shown in the 
plot. Additionally, the relation for early and late-type galaxies taken 
individually are independent on local galaxy density, implying that 
the increasing of M+/L r as the environment becomes denser is re- 
lated to the prevalence of massive systems in high-density regions, 
as shown in Fig. [5] Only for the most luminous galaxies there is an 
increment of M+/L r as a function of local density. 



5.5 Principal component analysis 

In the previous sections we have investigated the environmental de- 
pendence of the mean light-weighted stellar age, stellar mass, stel- 
lar metallicity and mass-to-light ratio in order to understand how 
galaxy properties are related to the environment. Here we will in- 
vestigate this issue through a different approach, based on an useful 
statistical technique frequently adopted to finding patterns in high- 
dimensional data. 

The parameter set discussed in this section was examined in 
more detail by means of a principal component analysis (PCA) 
in order to determine the combinations of parameters which sum- 



marise the distribution of the whole set. PCA is able to identify 
patterns in a data set and express the data in such a way as to 
highlight their similarities and differences. This technique is par- 
ticularly useful for reducing the dimensionality of a data space by 
identifying the linear combinations of input para meters with max- 
imum variance, called 'prin cipal components' l IMurtagh & Heckl 
ll987l ; ISodre & Cuevas|[l997h . 

A PCA was carried out on the combination of the parame- 
ter set { (log i* ) l , log Mi, , log Eio , log L r , log M+jLr } to deter- 
mine the variables that account for the large amount of variance 
and hence provide a useful description of the environmental depen- 
dence of galaxy properties. We have followed the same procedure 
as used in the previous sections by dividing the galaxy sample in 
bins of luminosity and stellar mass containing the same number of 
objects. A PCA was then executed for each bin separately. 

We found that the principal component of highest variance 
(PCI) accounts for more than 90 per cent of the total variance for 
all bins considered. The remaining variance is mainly accounted 
for the second principal component (PC2). Although each princi- 
pal component is a linear combination of all the input variables, 
it is interesting to verify whether some variables in the parameter 
set are well correlated with the main components because, in this 
case, they are responsible for a significant fraction of the sample 
variance. In Table|2]we show the Spearman rank correlation coeffi- 
cients between the two principal components and the five parame- 
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ters in our input set, for each bin divided by luminosity and stellar 
mass. We found that the local galaxy density (log Eio) shows the 
strongest relation with PC 1 ; the Spearman rank correlation coeffi- 
cient between this parameter and PCI is almost unity. Moreover, 
the PC2 shows a significant correlation with the {logt*)^, with 
correlation coefficient greater than 0.9. Thus, the plane PCI x PC2 
is well represented by the age-density relation in all bins of lumi- 
nosity and stellar mass investigated here. 

Additionally, to verify the behaviour of the other parameters 
we are investigating, we have excluded the mean stellar age from 
the input set. The PCA was then carried out on the parameter set 
{log Ah, log Eio, log L r , M+/L r }. For this set, the results indi- 
cate that the first principal component alone explains more than 
95 per cent of the sample variance for bins divided by luminosity 
and stellar mass. The second principal component is reponsible for 
almost all the remaining variance. The Spearman rank correlation 
coefficients among these two principal components and the param- 
eters in the input set are shown in Table[3]for each bin divided by lu- 
minosity and stellar mass. We found that the local galaxy density is 
again well correlated with the first principal component (rs = 1). 
The second principal component correlates well with both M* and 
M+fLr for bins in L r and with log L r and M±/L r for bins in M*. 

These results give clear evidence of the SFR-density relation 
present in our data since the local galaxy density is the primary 
responsible for the largest amount of sample variance and the mean 
light- weighted stellar age of galaxies is the secondary source of that 
variance. Moreover, we note that galaxy luminosity does not have 



any significant correlation with the principal components, whereas 
stellar mass and mass-to-light ratio can be considered as the drivers 
of the sample variance if one excludes the SFR-density relation. 



6 DISCUSSION 
6.1 Overview 

We have investigated the environmental dependence of some stellar 
population properties of galaxies in order to advance in our current 
understanding about the evolution of galaxies in distinct environ- 
ments. The most intriguing result shown here is that the relations 
between local galaxy density and properties related to star forma- 
tion is distinct when we split galaxies in bins of luminosity or stellar 
mass. 

The results shown in the last section will be inspected here in 
the light of two evolutionary paths by which galaxies can evolve. 
In the first one, galaxy properties (mainly related to star forma- 
tion and gas properties) are affected by environment through well 
known physical mechanisms acting on galaxies. This path, linked 
directly to the environment, gives origin to a nurture perspective 
for galaxy evolution. The second path is related to the initial condi- 
tions established during galaxy formation, which could account for 
the relations between galaxy properties and environment. Thus, it 
is related to a nature perspective driving galaxy evolution. 
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Figure 12. Median values of the mass-to-light ratio as a function of local galaxy density for galaxies in bins of luminosity. The plot in (a) is for all galaxies in 
our sample, (b) for early-type galaxies with D ?l (4000) > 1.67 and (c) for late-type galaxies with D ?l (4000) < 1.67. The numbers at the top-right are the 
count of galaxies in each panel. 
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Table 2. Spearman rank correlation coefficients for the five parameters set. 



Table 3. Spearman rank correlation coefficients for the parameter ser ex- 
cluding the mean stellar age. 



6.2 The star-formation-density relation: a nurture 
perspective? 

Recent galaxy redshift surveys have improved our understanding 
about what is happening with the star formation in galaxies in- 
habitin g different environments. Work s carried out b y Lewis et al.l 
J2002h using the 2dFGRS data, and iGomez et all J2003T) using 
the EDR SDSS data, have revealed a characteristic density at 
~ 1 ^75 Mpc~ 2 (corresponding to a cluster-centric radius of 3-4 
virial radii) associated to a "break" in the SFR-density relation. Be- 
low this density the SFR increas es only slight l y, whe reas at denser 
regions it is strongly suppressed. iTanaka et al.l J2004T) have comple- 
mented these studies by finding that the break in the SFR-density 



relation is seen only for fainter galaxies (M* + 1 < M r < M*+2), 
while the relation for brighter galaxies (M r < M* + 1) shows no 
break. Various physical mechanisms are advocated to explain these 
trends, but in general it is supposed that star formation properties 
of bright galaxies are affected by low-velocity interactions, while 
those of faint galaxies are affected by starvation (or strangulation). 
In this work we have investigated relatively bright objects 
(M r ^ M* + 1.5) and for a set of galaxies with Ha line in emission 
we have found no break or characteristic density in the relation be- 
tween EW(Hq) and local galaxy density, independently of galaxy 
luminosity (see Fig.[6j. Apparently, this resu l t is in conflict with th e 
conclusions reached by iLewis et al.l g2002f) . IGomez et alj J2003I) . 
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and lTanaka et al.l J2004I) that the star formation properties of (faint) 
galaxies change abruptly when one goes to denser environments. 
Here we clarify this question by investigating the colour-density 
relation for galaxies in our sample . For comparison pu rposes, this 
procedure is similar to that done bv lTanaka et alj J2004I) and shown 
in their Fig. 2. It is worth stressing that galaxy colours are measured 
in an uniform way by the SDSS photometric pipeline, that is, they 
do not depend on the methodology employed to com pute them. In 
additio n, the colour-density relation, as discussed bv lTanaka et al.l 
( [20040 . is very similar to the EW(Ha)-density relation. 

We show in Fig.[l3]the relation between the (g — i) colour and 
the local galaxy density for galaxies divided in two bins of r-band 
luminosity. In the left panel, we show the relation for all galaxies 
in our sample. A clear break is seen at Eio ~ 1 /170 Mpc -2 in 
the relation for fai nt galaxies, in comp lete agreement with the re- 
sults discussed by iTanaka et al.l Q2004J). On the other hand, bright 
galaxies tend to become redder in dense environments following a 
monotonically increasing relation. In the right panel of Fig.|13|we 
show the same relation but now restricting the sample to a set of 
galaxies with Ha detected in emission (avoiding to include AGN 
hosts, as discussed in Section [3). In this case, galaxies also tend 
to become redder in denser environments, but without showing 
any break at a characteristic density, even for faint galaxies. Thus, 
from the inspection of Fig. Q~3] we conclude that the break seen 
in the col our-density relation , or in the SFR-density relation dis- 
cussed bv lGomez et alJJ2003h . is associated to the reduced fraction 
of star-forming galaxies in denser environm ents, as discussed i n 
Section [5~Tj confirming the resu lts shown b y Balogh et al.l 020041) . 
Thus, in lGomez et alj J2003h and lTanaka et al.l 1 12004ft . for instance, 
the inclusion of absorption-line dominated galaxies (those with 
negative EW values) originated the break in the star formation- 
density relation at a particular local density, which actually is re- 
lated to an environment where the fraction of galaxies without re- 
cent star formation activity begins to dominate. 

These results also indicate that the decrease in the fraction of 
star-forming galaxies in dense environments occu rs in a wide range 
of den sities. A similar conclusion was reached bv lMateus & Sodrei 
1 120041) in a study on the environmental dependence of the fraction 
of star-forming and passive galaxies in a sample of field galaxies 
drawn from the 2dFGRS data. Even in rarefied environments the 
fraction of star-forming galaxies decreases significantly with in- 
creasing local galaxy density. Thus, there is a general consensus 
that physical mechanisms suppressing star formation in galaxies 
are not inherent only to the densest environments related to galaxy 
clusters. It is also interesting to note that faint and bright galax- 
ies seem to have evolved through distinct evolutionary paths, and 
different physical mechanisms could be at work. 

However, the relevance of such mechanisms acti ng on galaxies 
today may be important only in some particular cases. iBalogh et al.l 
( 120041) have confirmed this trend by combining samples from the 
2dFGRS and SDSS, with emphasis on environments related to 
galaxy groups. These authors argue that galaxies in dense regions 
have had typically more interactions than those in low-density re- 
gions over a longer period of time. Thus, it seems that galaxy 
transformations induced by environ ment have taken place mo re ef- 
fectively at higher redshifts. In fact. iKauffmann et alj 020040 have 
demonstrated that relations between structural parameters and stel- 
lar mass have little dependence on local galaxy density, indicating 
that these relations were already established at early times. In this 
way, mechanisms that are effective in disturbing the structure of 
galaxies (affecting the concentration index and surface mass den- 
sity) are not favoured by these results. 



6.3 A natural path for galaxy evolution 

In the last section, we have shown that the nurture hypothesis which 
is invoked to explain galaxy evolution in distinct environments also 
explains the observed star formation-density relation in the local 
universe. Indeed, in the hierarchical scenario of galaxy formation, 
as galaxy clustering evolves the density around a galaxy tends to in- 
crease in all environments. Higher density probably means more in- 
teractions which could, in principle, play a significant role in orig- 
inating the relation between star formation activity in galaxies and 
local galaxy density. 

The first numerical simulation s with CDM initial condi- 
tions (e.g. iFrenk et al] Il985l Il988h already suggested that the 
morphology-density relation is a natural product of hierarchical 
clustering. The high-redshift progenitors of today cluster galaxies 
were formed in high-density regio ns, which tend to collapse ear- 
lier than low-density environments JKaiserll984l ; lDavis et aljl985l : 
iBardeen et al.lll98q) . Galaxies formed in high-density regions tend 
to be ellipticals, because disks are destroyed by frequent interac- 
tions and the hot gas does not cool enough to form new stars. Thus, 
it is expected that the oldest galaxies were formed earlier and pref- 
erentially in high-density regions; galaxies in low-density environ- 
ments are able to keep their cold gas and form stars. Consequently, 
a natural path for galaxy evolution emerges, with galaxies in denser 
environments being more evolved than those in low-density regions 
JBenson et all200ll). 

Recently. I Juneau et al.l ( 120051) . using a sample of high-redshift 
galaxies from the Gemini Deep Deep Survey (GDDS), investi- 
gated the star formation history of galaxies and its dependence on 
stellar mass. They showed that the most massive galaxies formed 
most of their stars earlier than intermediate and low-mass galax- 
ies. This trend gives support to an apparent anti-hierarchical be- 
havi our, commonly referred as 'downsizing ' in galaxy formation 
(e.g. ICowie et alj| 1996k iKodama et al.ll2004l) . This picture, where 
massive galaxies stopped star formation at early times, whereas 
most low mass systems continue to form stars acti vely, is supported 
by the analysis of galaxies in the local universe JKauffmann et al.l 
l2003al ; lHeavens et al.ll2004l ; ISEAGalli]) . as well as by the observa- 
tion of a large population of massive galaxies at high-redshift (e.g. 
iGlazebrook et alj|2004lchen & Marzkell2004 among others). 

These trend s have been recovered by recent high-resolution 
simulations (e.g. IWeinberg et al.ll2004f) and semi-analytic model s 
of galaxy formation (e.g. iMenci et alj|2005l : Ide Lucia et alj|2006l) . 



On the observational side, the discovery of an overdensity of galax- 
ies (or proto-cluster) at 2 = 2.3 by ISteidel et al.l ( 120051) . with 
a significant fraction of old galaxies, is consistent with the the- 
oretical expectation for the acceleration of structure formation in 
de nser regions. This is in agreement with the conclusions drawn 
bv lEinasto et al.l |2005), that clusters in high-density environments 
evolve more rapidly than those in low-density regions. 



6.4 A nature via nurture scenario 

In this work, we have confirmed that high-density environments 
are dominated by massive galaxies with the oldest stellar popula- 
tions in the local universe, which can be associated to the early-type 
galaxies of the morphology-density relation. We have also found 
that the recent star formation activity of less massive galaxies, as 
inferred by the mean light-weighted age of their stellar popula- 
tions, does not change significantly along the density range pre- 
sented in this work. This result suggests that these low-mass galax- 
ies have large amounts of star formation in low as well as in high- 
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Figure 13. Left panel: Median values (thick lines) and respective quartiles (thin lines) of (g — i) colour as a function of local galaxy density for all galaxies in 
our sample divided in two bins of galaxy luminosity. Right panel: The same, but now restricted only to galaxies with Ha detected in emission. The numbers 
at the top-right are the count of galaxies in each panel. 



density regions of the local universe. However, as the lowest stel- 
lar mass bin considered in this work is dominated by late-type 
galaxies, with the fraction of low-mass early-type galaxies being 
affected by the selection criteria adopted here, this result should be 
confirmed by other observations. At the other extreme, the most 
massive galaxies in our sample have high values of mean stel- 
lar age independently of the environment where they are found. 
In short, massive galaxies have older stellar populations every- 
where. Only for an intermediate mass range (~ 3 — 6 x 10 10 
Mq), a ssociated to the transition in galaxy properties (see, for in- 
stance, ISEAGal ij) . the recent star formation activity of galaxies 
appears to decrease with increasing galaxy density. On the other 
hand, when we take galaxies divided according to their luminosi- 
ties, instead of stellar masses, the star formation activity of brighter 
galaxies and of fainter ones decreases significantly when one goes 
from low to high-density environments, resulting in the well know 
star for mation-density relation for galaxies selected by luminos- 
ity (e.g. *~ 



Lewis et al 



ITanaka et a l. 2004; Rines et al. 2005) 



2002; Gomez et aT] |2003l ; iBalogh et alj|2004l ; 



We have also shown in Fig. [12] that, on average, galaxies 
with similar luminosities in high-density environments have higher 
mass-to-light ratios, that is, they tend to be more massive than 
their counterparts in regions of low density (see als o ITuIIvi 12005 



for a similar result related to galaxy groups, and lEinasto et al 



120051 for galaxy clusters). Thus, the fact that galaxies of all lu- 
minosities present star formation activity decreasing with increas- 
ing density is related to the higher stellar masses of galaxies in 
denser environments, independently of their luminosities. Addi- 
tionally, from recent results concerning the existence of a down- 
sizing in the processes regulating star formation in galaxies, we 
expect that most massive galaxies have stopped to form stars at 
early times, perhaps thanks to internally driven mechanisms such 
as a mass threshold below which star formation w as not success- 
ful JMartin & Kennicudl200l| ; I Jimenez et alj|2005l) . Thus, the star 
formation-density relation found for galaxies divided according to 
their luminosities, being steeper for fainter galaxies JTanaka et al.l 
I2004T) . should be a "natural" consequence of galaxy evolution. 

These results are also supported by the recent findings ob- 



tained by ITanaka et al.l J2005I) in a study on the build up of the 
colour-magnitude relation of galaxies as a function of the environ- 
ment. They also confirm the downsizing way for the evolution of 
galaxy properties, with the star formation processes being displaced 
from massive to low-mass galaxies, and from galaxies in denser 
to low-density environments, as the evolution proceeds. In addi- 
tion, our work can also be related to the recent study carried out by 
IPoggianti et al.l J20051) . who have investigated the evolution of the 
proportion of star-forming galaxies (defined with the EW [O II]) in 
clusters since z = 0.4 — 0.8 (from a sample obtained with the ESO 
Distant Cluster Survey) to 0.04 < z < 0.08 (from a SDSS cluster 
sample). Poggianti et al. argue that the star formation-density re- 
lation should have a 'primordial' component and an 'evolved' one, 
each of them showing distinct environmental dependences, and that 
the relation is established at very high redshift at the moment of 
formation of the first stars in galaxies. This scenario is particularly 
linked to the trends shown by our own results. 

The trend related to stellar metallicity shown in Fig. [9] rein- 
forces these ideas: low-mass galaxies in denser environments are 
metal-rich compared with those in low-density regions. This arises 
from the fact that denser regions tend to be more evolved than 
rarefied environments. Thus, low-mass galaxies in dense regions 
have already been established in such environments, but have not 
stopped to form stars until the present epoch, contrarily to the evo- 
lution of massive systems. In this sense, a natural path for galaxy 
evolution proceeds via a nurture way mainly at high-redshifts: mas- 
sive galaxies have been formed in denser regions and evolved in 
an accelerated way, contrasting with a more unsocial life of low- 
mass galaxies preferentially inhabiting low density regions of the 
universe. This is also related to the timing of the process of for- 
mation of early-type galaxies that occurred a bout 1-2 Gyr earlier 
in cluster environme nts compared to the field JThomas et al.ll2005l : 
IClemens et al.ll2006T) . In this view, denser environments behave as 
an accelerator of the galaxy formation process, amplifying small 
natural differences imprinted on galaxies in their formation epoch. 
This results that "nature" necessarily acts via "nurture" effects. 

In this nature via nurture scenario for galaxy evolution, the 
environmental dependence of star formation properties of galax- 
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ies, commonly related to some (or various) physical mechanisms 
acting on late-type galaxies infalling onto galaxy clusters, is de- 
termined by the higher M±jL ratio of galaxies in denser envi- 
ronments, in addition to the observed downsizing in galaxy for- 
mation. These two ingredients play a crucial role in defining the 
star formation-density relation observed in luminosity divided sam- 
ples. Nevertheless, we note that environmental effects on galaxy 
properties really take place in some case s, such as the existence 
of post-starburst galaxies in clusters (e.g. IDressle r & Gunn 1 19831 ; 
iDressler et all 1 19991: IPoggianti et al.|[l999l : Ifran et al.l 12003) . pas- 



JBalogh 



ct 



afl l 19991 : 



sive spirals (e.g. Goto et al. 20031), and s hort starburst galaxies 



iMateus & So dre 2004). However, the fraction 



of such 'environment-product' objects is very reduced, indicative 
of a very rapid and efficient mechanism acting on these galaxies 
or actually due to the rarity of such transformations in the nearby 



7 SUMMARY 

In this fourth paper of the SEAGal collaboration, we have shed 
some light on the discussion about the environmental dependence 
of galaxy properties in the local universe. We based our analysis 
on the stellar population properties of galaxies in a volume-limited 
sample drawn from the SDSS Data Release 4. The application of a 
spectral synthesis method to the data produces robust estimators for 
mean stellar ages, mean stellar metallicities and stellar mass, which 
have been used in this work to characterise the stellar populations 
of galaxies. The environment is defined by the local galaxy density 
estimated from a nearest neighbour approach, and by the distance 
to cluster centres obtained from a public catalogue of SDSS clus- 
ters. The main approach used in this study is the comparison of the 
relations between galaxy properties and environment for galaxies 
divided in intervals of luminosity and stellar mass. We summarise 
our main findings below: 

(i) We recover the star formation-density relation in terms of the 
mean light- weighted stellar age, (logi»)i, which is strongly cor- 
related with star formation parameters derived from Ha emission 
line. 

(ii) We confirm that high-density environments are populated by 
a large fraction of massive galaxies with old stellar populations, in 
opposition to low-density regions, dominated by low-mass galaxies 
actively forming stars. We also note that galaxies with intermediate 
stellar masses have constant fractions along the range of densities 
covered by our sample. The transition in galaxy fraction occurs at 
E* ~ 0.7 h^ Mpc -2 , corresponding to 2-3 virial radii from clus- 
ter centres. 

(iii) The environmental dependence of (log£*)i is distinct 
when we divide galaxies according to luminosity or stellar mass. 
The relation between mean age and local density is remarkable for 
galaxies in all bins of luminosity. On the other hand, only for a in- 
termediate stellar mass interval (associated to a transition in galaxy 
properties) the relation shows a change in galaxy properties. For 
low-mass galaxies, the relation slightly increases along the density 
range, but showing low values of (log t*)i,- characteristic of late- 
type galaxies - in all environments. The most massive galaxies also 
show a slight change in the median values of (logt*)x as galaxy 
density increases, but their values are higher everywhere. 

(iv) We find that the distinct behaviours of the age-density re- 
lation for galaxies divided by luminosity or stellar mass are as- 
sociated to the large stellar mass of galaxies with the same lumi- 
nosity in dense environments. In other words, the well known star 



formation-density relation results from the prevalence of massive 
systems in high-density environments, independently of galaxy lu- 
minosity, in addition to the observed downsizing in galaxy forma- 
tion, where the star formation is shifted from massive galaxies at 
early times to low-mass ones as the universe evolves. 

(v) A principal component analysis of our parameter set reveals 
that the local galaxy density is the primary responsible for the to- 
tal variance present in our data, whereas the mean light-weighted 
stellar age of galaxies is the second one. This result reflects that the 
age-density relation is the main driver of the environmental depen- 
dence observed in some galaxy properties. 

(vi) The mean stellar metallicity of less massive/luminous 
galaxies increases in high-density regions, indicating that even low- 
mass galaxies in dense environments tend to be more evolved than 
their counterparts in low-density regions. 

(vii) Our results support that a natural path for galaxy evolu- 
tion proceeds via a nurture way mainly at high-redshifts: massive 
galaxies have been formed in denser regions and evolved in an ac- 
celerated way, contrasting with a more unsocial life of low-mass 
galaxies preferentially inhabiting low density regions of the uni- 
verse. 



In this work, the concept of 'nature via nurture' (see iRidlevI 

120031 for a genetic view of this expression) in the scenario of galaxy 
evolution was advocated in order to summarise the results obtained 
here related to the environmental dependence of galaxy properties 
in the local universe. Many steps are needed towards a compre- 
hensive view of the processes which have guided galaxy formation 
and evolution, mainly those related to the properties of the stellar 
content of galaxies in high redshifts. We expect to contribute fur- 
ther to these issues in other papers of this series on Semi-Empirical 
Analysis of Galaxies. 
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